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Abstract: The Great Xing’an Mountains boreal forests were focused on 
in the northeastern China. The simulated future climate scenarios of 
IPCC SRES A2a and B2a for both the baseline period of 1961-1990 and 
the future scenario periods were downscaled by the Delta Method and the 
Weather Generator to produce daily weather data. After the verification 
with local weather and fire data, the Canadian Forest Fire Weather Index 
System was used to assess the forest fire weather situation under climate 
change in the study region. An increasing trend of fire weather severity 
was found over the 21st century in the study region under the both future 
climate change scenarios, compared to the 1961-1990 baseline period. 
The annual mean/maximum fire weather index was predicted to rise 
continuously during 2010-2099, and by the end of the 21st century it is 
predicted to rise by 22%-52% across much of China’s boreal forest. The 
significant increases were predicted in the spring from of April to June 
and in the summer from July to August. In the summer, the fire weather 
index was predicted to be higher than the current index by as much as 
148% by the end of the 21st century. Under the scenarios of SRES A2a 
and B2a, both the chance of extremely high fire danger occurrence and 
the number of days of extremely high fire danger occurrence was pre¬ 
dicted to increase in the study region. It is anticipated that the number of 
extremely high fire danger days would increase from 44 days in 1980s to 
53-75 days by the end of the 21st century. 
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Introduction 

As one of the largest terrestrial biomes in the world, the boreal 
forest, which accounts for approximately 30% of the world’s 
forested areas, plays an important role in the global cycling of 
energy, carbon and water (Conard and Ivanova 1997). The study 
region in the Great Xing’an Mountain Range (GXAMR) is the 
most important boreal forest in China, and it is also the most 
southern part of the global boreal forest biome. 

Forest fires are strongly influenced by weather and climate 
(Flannigan and Harrington 1988; Johnson 1992; Swetnam 1993). 
Therefore, forest fire is one of the most important climate change 
impacts on the boreal forest now and in future because of its 
dominant role in interacting with other impacts (Wheaton 2001). 
With approximately 0.15 million hectares burned annually, the 
fire occurrence in GXAMR is relatively less extensive than the 
average level of the global boreal forest (Flannigan et al. 2005). 
Due to the proactive fire suppression policy carried out by the 
forest administration in the region for the past two decades, all 
forest fires were put out immediately after the ignition was de¬ 
tected. As a result, fuel was rapidly built up in the study region, 
and it became a serious concern as fuel means with higher fire 
risk and more severe damage. The accumulation of deadwood 
and debris across the continuous boreal forest landscape in 
GXAMR could result in severe forest fires in the region, which 
may be too difficult and too costly to suppress (Andrews et al. 
2007). 

It is generally accepted that the global temperature is increas¬ 
ing with the largest increase happening at high latitudes in the 
Northern Hemisphere where boreal forest resides (Soja et al. 
2007). Climate change projection indicates that GXAMR is 
likely to become hotter and drier in the future; therefore, un- 
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der-standing interactions between the climatic influence and the 
characteristics of forest fire regimes is essential for fire man¬ 
agement in the region. Studies demonstrates that forest fires are 
influenced by the changing climate (Westerling et al 2006; Ly 
and Tian 2007; Lavorel et al. 2007; Running 2006; Pausas 2004; 
Mollicone et al. 2006). A number of studies focused on the pos¬ 
sible impacts of global warming on wildland fires using the 
General Circulation Models (GCMs), which can be used to pro¬ 
ject future climate change, the impact of this change on fire 
weather severity (Flannigan et al. 2005; Pino et al. 1998; Stock et 
al. 1998; Hennessy et al. 2006), and potential changes in areas 
burned due to climate change (Stocks et al. 1998; Tymstra et al. 
2007; Flannigan et al. 2005). All these studies were conducted in 
temperate regions of the western US, Australia and the boreal 
forests in Canada or Russia. As an important component of the 
global boreal forest, the boreal forest in China’s GXAMR region 
needs a similar study to assess possible associations between 
forest fire and climate change. 

The Canadian Forest Fire Weather Index System (CFFWIS) is 
used to quantitatively measure fire risk for the study region. With 
daily temperature, relative humidity, precipitation, and wind 
speed as input, CFFWIS can be used to keep track of the forest’s 
day-to-day susceptibility to fire (Van Wagner 1987). The system 
is comprised of several standard output components. The first 
three are fuel moisture codes that follow daily changes in the 
moisture contents of three layers of forest fuel with different 
drying rates. The last three components are the fire behavior 
indexes representing the rate of fire spread, fuel weight con¬ 
sumption and fire intensity. The system depends solely on 
weather readings taken each day on temperature, relative humid¬ 
ity, wind speed and precipitation (if any) during the previous 24 
hours. Three fuel moisture codes are the Fine Fuel Moisture 
Code (FFMC), Duff Moisture Code (DMC), and Drought Code 
(DC). The Fire Weather Index (FWI) is derived from the fuel 
moisture codes and weather readings. FFMC, DMC, DC and 
FWI are arranged in code form with values rising as the moisture 
content decreases or the fire danger increases. There are two 
other outputs from the system: the Daily Severity rating (DSR) 
and the Seasonal Severity rating (SSR). In general, FWI reflects 
the fire danger and is usually used as fire hazard warning. DSR is 
a measure on the difficulty of fire control (Di 1993), and SSR is 
an index representing the seasonal mean of fire control difficulty 
(Flannigan et al. 2000). FWI is adopted as a comprehensive 
measure of forest fire danger in the study region. 

The purpose of this study is to assess potential changes in for¬ 
est fire weather risk influenced by climate. This change in risk is 
related to forest fire activities such as fire occurrence, area 
burned, and the impact of such changes on fire activities in the 
GXAMR region. Such information is very critical for medium to 
make long term forest management plan in the region. To 
achieve this goal, the weather data collected from the region’s 
weather stations from 1950s to 2000s were used. Then, GCMs 
are used to simulate future weather conditions until the end of 
21st century under two climatic change scenarios of the Hadley 
Center Coupled Model 3-A2a and B2a (Gordon et al. 2000; 
Pope et al. 2000). The simulated daily weather conditions such as 
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daily temperature, relative humidity, wind speed and precipita¬ 
tion are then used as input to CFFWIS to generate FWI, which 
reflects daily fire risk. The simulated FWI provides a long term 
trend of forest fire potentials in the region, and is an important 
reference for assessing the region’s forest fire risk under the 
global climatic change scenarios, most likely to occur in the 
future. 

Materials and methods 

Study area descriptions 

The GXAMR region is located in the northeastern China 
bounded by 50°10'N-53°33'N and 121°12'E-127°00'E with a 
total area of 8.35 million ha (Fig. 1). The region locates in the 
cool temperate and a terrestrial monsoon climate zone with mean 
annual temperature of -2.8°C and mean annual precipitation of 
746 mm. The forest landscape in this region is dominated by 
Xing’an larch ( Larix gmelini), and other major tree species are 
Asian white birch ( Betula platyphylla), Scot pine ( Pinus sylves- 
tris var. mongolica) and Mongolian oak ( Quercus mongolica ) 
(Heilongjiang Yearbook Editorial Board 2008; Jiang et al. 2002). 

Data collection 

Meteorological data 

Meteorological data are collected from 18 weather stations in 
and around the study region, and the data comprises of daily, 
monthly and annual observations on temperature, wind speed, 
relative humidity and precipitation. Among 18 weather stations, 
seven are located inside the study region, and the other 11 sta¬ 
tions are in the surrounding area (Fig. 2). Interpolation is used to 
generate grid data of weather factors for the later modeling exer¬ 
cise. The purpose of selecting the data from the surrounding 
areas is to make the interpolated results more accurate. The 
weather data was downloaded from China’s meteorological data 
sharing service system. The availability of the weather data is 
listed in Table 1, in which there are the time and type of the data 
available. The data are not complete for all 18 weather stations, 
and some data are missing in several stations for certain time 
periods. For example, there are missing data in Tahe station for 
1962-1971 and in Mohe station in 1958. Therefore, 
one-dimensional linear regression method has to be used to fill in 
the data gaps. 

Climatic scenarios 

The data for the climate change scenario in this study include the 
observed mean monthly values of mean/maximum temperature, 
total precipitation, wind speed and relative humidity for 
1961-1990 (referred to as 1980s). The projected mean monthly 
values are used for 2010-2039 (referred to as 2020s), 2040-2069 
(referred to as 2050s), and 2070-2099 (referred to as 2080s). The 
four periods have the same length of 30 years. The projection is 
done by using the third version of the Hadley Center Coupled 
Model (HadCM3). As the most commonly used projection 
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method, HadCM3 requires no flux adjustment or the stable cli¬ 
mate in global mean, and it demonstrates to perform better than 
other models for Eastern Asia, where China is located. The sce¬ 
narios tested belong to the A2 and B2 storylines and the families 
of A2a and B2a, respectively. In particular, the A2a storyline 


family describes a very heterogeneous world, presuming 
self-reliance and preservation of local identities. The B2a story¬ 
line emphasizes the local solutions to economic, social, and en¬ 
vironmental sustainability (Leckebusch and Ulbrich 2004). 



MONGOLIA 


Great Xing'an 
mountains 


Q Boreal Forest 

International Boundarl 


e Circumpolar Extent of the Boreal Forest 


Inner Mongolia 


Heilongjiang 


Fig. 1 (a): the geographic location of the study region; (b): the northern circumpolar of the global boreal forest. 
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Fig. 2 The grid boundary for conducting the HadCM3 model 
analysis that covers the study region and 18 weather stations. The 

grid of interpolation has 1000 rows and 1000 columns with the spatial 
resolution of 40 m x 40 m. 

Statistical analyses 

Delta change methods 

Generally, the spatial resolution of the output from GCMs like 
HadCM3 is coarse and not suitable for representing regional 
climate/weather situations. Therefore, downscaling is needed to 
make the resolution finer to meet the regional scale requirement. 


In the last decade, downscaling techniques for both dynamic and 
statistical approaches were developed. The coarse spatial resolu¬ 
tion can thus be made finer in the studies of climate change sce¬ 
narios. 


Table 1. Data type and availability 



The data type and the time period when the data available 

Weather - 

variable 

Annual, 

1951-2008 

Monthly, 

1961-1990 

Daily, 

1951-2008 


Mean T 

Mean T 

Maximum T 

T(°C) 

Maximum T 

T anomaly 

Maximum T 


RH (%) 

Mean RH 

Mean RH 

Mean RH 

P (mm) 

P 

P anomaly (%) 

P 

Number of days 

with P > 0.1 mm 

20-20h P 

W (m/s) 

Mean W 

Mean W 

Mean W 


T: Temperature; RH: Relative humidity; P: Precipitation; W: Wind 
speed. 

The Delta change method (Akhtar et al. 2008; Zhao and Xu 
2007) is recommended by the Global Change Research Program 
of the United States. This is a simple and useful method to gen¬ 
erate future climate variables by combining previous observa¬ 
tions and the simulated results through the HadCM3 model. The 
method is successfully used in many studies on the impact of 
climate change (Akhtar et al. 2008). The weather data used as the 
input to the Delta change method are the observations of tem¬ 
perature, relative humidity, wind speed and precipitation from 18 
weather stations covering from 1961 to 1990. The future monthly 
values of the four weather variables were constructed using Eq. 
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1-4 as follows for each grid cell in the study region: 

T, = T + (f, - f ) (l) 

p,=pjf,p;') < 2 ) 

r,=r 0 {r,K') < 3 > 

w, =w,(w,w;') < 4 ) 

Where, Tf is the monthly temperature for 12 x 30 months in each 
of the three future scenario periods (2020s, 2050s and 2080s). T 0 
is the observed monthly temperature for 12 x 30 months in the 
baseline period of 1961-1990. 7T is the mean of the future 

monthly temperature generated from the HadCM3 model for the 
future scenario periods, and T o is the mean of the monthly tem¬ 
perature for the baseline period. The other three equations are 
used for deriving monthly precipitation, relative humidity and 
wind speed for the future scenario periods. The variables in the 
equations with the same subscripts as the ones in Eq. 1 denote 
the same time period. 

Weather generator 

CFFWIS requires the daily weather conditions of temperature, 
precipitation, relative humidity and wind speed as input, and the 
weather data coming from the Delta method is monthly. Hence, 


we need to generate daily weather data based on the monthly 
information for the future scenario periods. Weather generators 
are computational tools to produce daily and site-specific climate 
change scenarios (Semenov et al. 1998). One of the most com¬ 
monly used weather generators is WGEN by Richardson and 
Wright (1984). It is chosen to produce daily weather data for 
each of the future scenario periods. WGEN produces synthetic 
daily time series of a group of climate variables such as maxi¬ 
mum temperature, precipitation, relative humidity and wind 
speed with certain statistical properties. The monthly temperature, 
precipitation, relative humidity and wind speed are among the 
input to WGEN. The observed daily data and the daily data gen¬ 
erated through WGEN for the baseline period were compared 
and tested at each weather station, and the test results showed 
there are low errors between the two data sources. Also, the se¬ 
rial correlation and cross serial correlations among the generated 
data are acceptable. For example, the relationships between the 
observed and the WGEN generated temperature and relative 
humidity for the baseline period for the Mohe weather station are 
shown in Fig. 3, showing a strong correlation between the ob¬ 
served and generated data. To convert point data of 18 weather 
stations to the grid data covering the study region, the ordinary 
Kriging analyses of the Geostatistical Analyst extension of the 
ESRI ArcGIS software (Yost et al. 1982; Chien et al. 1997) is 
applied in the study, and it is the most flexible Kriging model 
because it assumes a constant but unknown mean. 



Observed temperature (°C) Observed relative humidity (%) 


Fig. 3 (a): Scatter plot of observed and simulated daily maximum temperatures for the Mohe Weather Station, (b): Scatter plot of observed 

and simulated daily relative humidity for the Mohe Weather Station. The simulated data are based on the HadCM3 model and the WGEN model and 


for 1966-2008. 


Testing FWI 

FWI was used to link weather and climate conditions to fire haz¬ 
ards (Wheaton 2001). The system was once used in the study 
region during 1981-1991, but other studies suggested the system 
could be directly applied to boreal forests without the need for 
local calibration (Stocks et al. 1998; Willis et al. 2001; Tian et al. 
2006). FWI is used to make daily fire hazard projections for each 
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of the three future scenario periods based on the daily weather 
information generated from the WGEN method. Before FWI is 
used, the system is first tested by using both local weather data 
and fire data, which are provided by the local forest administra¬ 
tion. The fire data for testing covers that in 1966-2008, including 
fire ignition time, ignition location, cause and area burned. The 
test results indicate that a significant correlation (r = 0.91) exists 
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between the annual mean FWI and annual fire occurrences (Fig. burned is poor. 
4), although the correlation between the FWI and annual area 



Fig. 4 (a): Scatter plot of the mean daily FWI and the mean daily fire occurrences versus Julian day for 1966-2008. (b): Scatter plot of the 

mean daily fire occurrence versus the mean daily FWI values for 1966-2008. 


Results and discussion 

In the study region, annual mean temperature increased by ap¬ 
proximately 2.0°C since the 1950s (Fig. 5a), which is substan¬ 
tially greater than the global average increase. HadCM3 reflects 
that the annual mean temperature will potentially increase by 
4.7°C under the A2a scenario, and it also reflects that the annual 
mean temperature will increase by 3.9°C under the B2a scenario 
during 2010-2099 in the study region. The projected precipita¬ 
tion also showed an increasing trend during the period (Fig. 5b): 
a 30% increase was under the A2a scenario, and a 16% increase 
was under the B2a scenario. It was observed that since 1959 the 
annual mean temperature anomaly exhibited an increasing trend 
during the fire season (April-October), with the highest positive 
anomaly occurring by the end of the 1980s (Table 2). Coupled 
with the rising temperature anomaly, the annual precipitation 
anomaly percentage also shows an increasing trend during the 
fire season for 1959-2008, although the increase is slight (Table 
2 ). 


Table 2. Summary of the proportions of fire activities corresponding 
to the fire danger classes during 1966-2008. 


Time period 

TT 7 _ j 1 * 1 1 1 

weamer vanaoies 

I 

II 

III 

IV 

V 

Annual mean temperature 

anomaly (°C) 

-4.53 

-4.03 

-2.55 

4.72 

5.85 

Precipitation percentage 

anomaly (%) 

0.15 

-5.65 

-0.93 

6.95 

0.44 


I: 1959-1968; II: 1969-1978; III: 1979-1988; IV: 1989-1998; V: 
1999-2008 


In the study region, there were 36 fire occurrences burning 
0.15 million ha on average annually during 1966-2008. A large 
forest fire occurred in the region on May 6, 1987, which lasted 
for 25 days and burned 1.01 million ha, causing huge economic 
and ecological damage (National Bureau of Statistics and Minis¬ 
try of Civil Affairs 1995). When the fire data between the two 
periods of 1966-1987 and 1988-2008 were compared, both fire 
occurrence and burned area showed significant decrease after the 
1987 large fire incident. The average of annual fire occurrence 
during 1988-2008 decreased by 31.4% compared to that in 
1966-1987, the annual average burned area decreased by 77.7%, 
and the annual average burned forest decreased by 86.7%. How¬ 
ever, in 1988-2008, the overall fire activities became more se¬ 
vere, with considerable variation in the magnitude of both fire 
occurrence and area burned (Fig. 6a). This may be attributed to 
the warming climate. During the fire season, fire activities con¬ 
centrated in the spring months of March to May, and in the au¬ 
tumn months of September and October (Fig. 6b). However, the 
fire trend changed after entering the 21st century, with more fires 
occurring between July and August during 2001-2008. For ex¬ 
ample, 33.1% of all fires occurred during the summer months in 
this period, while only 3.8% of all fires occurred in summer in 
1966-2000. 

The emphasis of this study is on FWI and its changing trend 
under different climate change scenarios for the three future pe¬ 
riods. Although the projected FWI is a daily time series for the 
three future time periods, it is also examined at annual and 
monthly scales. The annual mean FWI was calculated at all 
weather stations and interpolated to generate maps as shown in 
Fig. 7. Under the scenario A2a, the annual mean FWI will in¬ 
crease by approximately 5% in the 2020s, 18% in the 2050s, and 
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46% in the 2080s, compared to that in the baseline period. Under 
the scenario of B2a, the annual mean FWI will decrease by 3% 
in the 2020s, but increase by 24% in the 2050s, and by 26% in 
the 2080s, compared to that in the baseline period. The annual 
maximum FWI was calculated, and the results indicate that it 
will continue to increase for 2010-2099. In the 2080s, it will 
increase by as much as 52% under Scenario A2a and 22% under 
Scenario B2a. As for the monthly time scale, almost all monthly 
mean FWI demonstrates an increasing pattern under the two 



climate change scenarios for the three future periods compared to 
that in the baseline period (Fig. 8). The plots in the figure show 
the peaks at May and June, but the highest increase of the FWI 
happens in July and August. For example, under Scenario A2a, 
the monthly mean FWI will increase by 29% to 41%, compared 
to that in the 1980s. The same FWI will increase to 148% by the 
end of the 21st century. April will be subject to the second high¬ 
est increase of the monthly average FWI, with the FWI going up 
by 50% in the 2080s. 



Fig. 5 (a): Observed (1959-2008) and projected (the three future periods under the two climate change scenarios) annual mean temperatures, 

(b): Observed (1959-2008) and projected (the three future periods under the two climate change scenarios) annual precipitations. 
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Fig. 6 (a): Time series of annual fire occurrence and area burned based on the observed data for 1966-2008. (b): Time series of monthly fire 

occurrence and area burned base on the observed for 1966-2008. 

It is customary to use an ordinal variable to rank fire danger on forming the general public. Fire danger was ranked into five 

different levels from low danger to extreme high danger for in- categories according to the FWI values as presented in Table 3. 
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The ranking is based on FWI values, and the fire activities cor¬ 
responding to the FWI ranks are also presented. Based on the 
ranking scheme, the daily FWI value is calculated for all three 
future periods under the two climate change scenarios, and it is 
classified into five fire danger classes (Fig. 9). In general, the 
number of low fire danger days will decrease and the number of 
extremely high danger days will increase for the three future 


periods of 2020s, 2050s and 2080s under both climate change 
scenarios. In the 1980s, on average there were 37 days of very 
high fire danger, and the number of days is predicted to double 
by the end of the 21 st century. In the entire baseline period, there 
were a total of 44 days of extremely high fire danger; and in the 
2080s this number will increase to 75 days under the A2a sce¬ 
nario and to 53 days under the B2a scenario. 
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Fig. 7 Maps of annual mean FWI values for the baseline and the 
three future periods under the two climate change scenarios. Map 
(a): baseline; Map (b): 2020s A2a; Map (c): 2050s A2a; Map (d): 
2080s A2a; Map (e): 2020s B2a; Map (f): 2050s B2a; Map (g): 
2080s B2a. 


Conclusion 

In the study region, the annual mean temperature increased by 
approximately 2.0°C in the past 50 years, and the increase is 
likely to intensify in the 21st century. According to the HadCM3 


model, the annual precipitation would increase by 16% to 30% in 
2080s under both climate change scenarios. Since 1959, the 
anomaly of annual mean temperature in the fire season exhibied 
a large increase, and the number of extremely hot days also 
raised. In contrast, the precipitation showed a slight increase as 
indicated by the anomaly of annual precipitation percentage. 
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Forest fires and their characteristics are strongly influenced by 
weather and climate. In the study region, it was noted that in¬ 
creasing climate warming substantially increased fire activities 
such as fire occurrence, area burned and fire size since 2000. 
Although most of the forest fires occurred in spring and autumn 
in 1966-2000, the recent evidence indicated that the fire pattern 
changed and more forest fires occurred in the summer since 2000: 
33.1% of fires occurred in summer since 2000 and 3.8% in sum¬ 
mer during 1966 and 2000. 
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Fig. 8 Mean monthly FWI values for the baseline and three future 
periods under the two climate change scenarios. 




Time period 


Fig. 9 Number of days with different fire danger ranking for the 
Springer 


baseline and three future periods under the two climate change sce¬ 
narios. 


Table 3. Summary of the proportions of fire activities corresponding 
to the fire danger classes during 1966-2008. 


Fire danger 

class 

FWI 

interval 

Proportion of fire activities (%) 

Occurrence 

AB 

FAB 

Low 

[0, 2) 

4.30 

0.66 

0.34 

Moderate 

[2, 5) 

11.30 

6.86 

1.91 

High 

[5,9) 

19.91 

11.51 

11.15 

Very high 

[9, 16) 

32.95 

38.27 

46.90 

Extremely high 

[16, ») 

31.54 

42.71 

39.71 


AB: Area burned; FAB: Forested area burned. 

In the study region, both annual mean FWI and annual maxi¬ 
mum FWI are predicted to increase by 22%-52% across much of 
the study region in 2080s under either of the two climate change 
scenarios. The monthly mean FWI suggests that the study region 
is going to experience the highest fire danger from March to 
November in the future periods rather than in the baseline period 
when the fire danger was highest in spring and autumn. The re¬ 
sults also indicate that the monthly mean FWI value could in¬ 
crease by as much as 148% by the end of the 21st century, and 
the increase would happen in spring and summer. Forest fires 
that occur on days with high fire danger tend to be large and 
severe fires. The study indicates that in the future the number of 
days with extremely high fire danger is going to increase (53-75 
days of extreme high danger by 2080s, compared to a present 44 
days of the same danger rating). As a result, the proportion of 
large fires among all fire sizes in the future is very likely to in¬ 
crease. 

Under either of climate change scenarios, the FWI in China’s 
boreal forest is going to increase in the three future periods, 
which indicates that the forest fire danger is going to increase. 
Such information is an alert to the forest management agencies in 
the study region. The information on other factors influencing 
fire ignition and fire behavior such as terrain, vegetation, fuel 
type and load, topography, transportation network, fire suppres¬ 
sion capability, and so forth, should also be helpful to forest 
management to set up their protection plan. At present the den¬ 
sity of meteorological stations is relatively low, the spatial dis¬ 
tribution of the stations is uneven and some of the stations are 
not in the forested area. These factors affect the accuracy of the 
FWI derivation. Nevertheless, the study results provide impor¬ 
tant information on the relationship between future fire activities 
and the impact of climate change on forest fires in China’s boreal 
forest. Our results could also be a helpful reference for similar 
studies in other boreal forests around the world. 
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